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Abstract 

We investigate the energy landscape of two dimensional network models for 
covalent glasses by means of the lid algorithm. For three different particle 
densities and for a range of network sizes, we exhaustively analyse many 
configuration space regions enclosing deep-lying energy minima. We extract 
the local densities of states and of minima, and the number of states and 
minima accessible below a certain energy barrier, the 'lid'. These quantities 
show on average a close to exponential growth as a function of their respective 
arguments. We calculate the configurational entropy for these pockets of 
states and find that the excess specific heat exhibits a peak at a critical 
temperature associated with the exponential growth in the local density of 
states, a feature of the specific heat also observed in real glasses at the glass 
transition. 

61.43.-j,61.43.Fs,5.40.+j,5.70.Fh 
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I. INTRODUCTION 



Since the beginning of the last decade, systems with complex multi-minima energy land- 
scapes have attracted increasing attention with a common theme being thermal relax- 
ation or more generally, stochastic dynamics on the landscapes. Such dynamics can either 
have intrinsic physical interest or be utilized as an optimization device as done in annealing 
techniques. A number of approaches have been developed, focusing on different aspects 
of the problem: On the one hand, molecular dynamics and Monte Carlo simulations are 
performed, often using highly refined model potentials, which are designed to reproduce as 
closely as possible the actual dynamics at short times and the equilibrium statistical 
mechanical properties of the system []3|-[5| , respectively. On the other hand, one uses simple 
models describing only selected features of the system, which are amenable to analytical 
techniques ||[7| or can be studied numerically || in enough detail to yield general insights 
into the qualitative and semi-quantitative behaviour of the system. As part of the latter 
approach one can consider abstract graph models, which formally can be thought of as 
"lumped" representations of the energy landscape itself P,|9|JT0|. 

The network models for covalent glasses presented in this paper belong to the second 
class of approaches, since they are tailored to describe the slow part of the complex hierarchy 
of relaxational degrees of freedom JE^-H] which characterize glasses: On the shortest time 
scales, we are dealing with small vibrations in the immediate vicinity of individual minima 
of the energy surface. These are responsible for most of the vibrational and reversible elastic 
properties of the glass. Here, the analysis usually employs matrix diagonalisations at the 
point of the minimum, or short time MD- simulations. At the next level, neighbouring min- 
ima are accessible by crossing very small barriers. This mechanism is probably responsible 
for some of the anomalous low-temperature properties of glasses. One would suspect that, 
at this level of detail, the so-called "two-state- models" and their descendants e.g. the 



soft potential models |17|,18|, would be a relevant theoretical description, which can be com- 
plemented by studying the diffusion of (single) particles by MD/MC-simulations at various 
temperatures. At large time scales and/or high temperatures (up to the point where the 
glass melts), the main structural feature of the glass is its topology. Accordingly, one usually 
visualises the glass as a random network of building units [|1^,|14|,|]1| , where the links repre- 
sent either chemical bonds (e.g. Si — Si) or sequences of chemical bonds (e.g. Si — O — Si 
or B — O — B). The relevant excitations are likely to be long wavelength distortions of the 
covalent network, which involve the displacement of many atoms, with each displacement 
small compared to the interparticle distance. Such distortions can substantially change the 
geometry of the structure, while they only weakly affect its topology. In systems containing 
thousands of atoms per simulation cell it becomes computationally very expensive if not 
impossible to run MD/MC-simulations for the required times while still using highly refined 
potentials. However, the structural and energetic hierarchy of a covalent glass, leading to 
a separation into vibrational, geometric, and topologic properties of the glass, opens the 
possibility of employing network models on lattices to selectively describe the topology of 
the glass. Similar lattice models for polymers have been successfully analysed in recent years 
using MC-simulations [f|. 

A salient feature of glasses is the glass transition, with a peak of the (excess) specific heat 
capacity C p at a temperature near the transition temperature T G |T3|JT^| . This peak is usually 
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associated with the so-called configurational entropy reflecting the multitude of different 
topological structures accessible to the glass during this transition. Each configuration 
represents a basin around a relatively stable local minimum of the potential energy. Thus, 
the configurational entropy is an excess entropy of the glassy state relative to the crystalline 
state, whose entropy at this temperature is dominated by the vibrational states. Based on 
the previous discussion, it should therefore be possible to link the excess entropy of glasses 
to statistical features of the energy landscape of covalent network models. 

Motivated by these considerations we numerically analyzed the energy landscape of small 
network models on two-dimensional lattices for a range of sizes and densities. The restriction 
of the nodes of the networks to a lattice allows complete characterization of subsets of 
the discrete landscape. Using the lid algorithm ||, we performed exhaustive searches of 
local regions around deep-lying minima (so-called pockets), which yield information on the 
local densities of states, the available configuration space volume and the distribution of 
neighbouring minima. This information is then used to understand some of the features of 
the thermodynamics and dynamics of the system. 



II. MODEL AND ALGORITHM 

A. Lattice network model 

The networks were placed on square lattices with periodic boundary conditions. The size 
of the repeated cell S x S ranged from 10 x 10 to 20 x 20 grid points. The number density 
p = Na/V, V = S 2 , of the Na building units per cell ( such units are henceforth for simplicity 
called atoms) was chosen to be approximately 0.13, 0.14 and 0.15. The interaction potential 
given by Eqs. and 0) between the atoms consisted of a sum of a two body and a three body 
term. The former, V^r), grows quickly towards large positive values for distances r < 1.8a 
and equals infinity for r < a, while for r > 3.2a it smoothly approaches zero. The lattice 
parameter a was chosen in such a way that the optimal distance between atoms was about 
2.1a. The three body term V^(9) details the angular dependence of the interactions among 
nearest neighbors (it only applies for r < 3.2a). It has a minimum at about 120°, reaches 
infinity for angles smaller than 80° , and vanishes smoothly when the angle approaches 
180°. The actual formulae used are not very important, but are nevertheless mentioned for 
completeness: 

{+oo : r < 1.6 
l/3(r-1.6)(r-3.2) 2 (r-7.1) : 1.6 < r < 3.2 (la) (1) 
: r > 3.2 

and 

9 < 80° 

80° < 9 < 180° (16) (2) 
9 > 180° 

Thus, four-fold coordination was possible, but not favoured energetically. A typical 
metastable configuration is shown in fig. 1. The binding energy of such deep-lying minima 



V 3 (9) = 5.8 x 10- 8 (fl - 80°) (9 - 180°) 2 9 
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was in the range of —6.0 to —6.5 eV/atom. The (crystalline) ground state of an infinite 
system - without an underlying lattice - has a hexagonal structure similar to that of a 
graphite layer, with each atom surrounded by three equidistant neighbours. In the finite 
on-lattice system, most of the deep-lying minima encountered contained a certain amount 
of adjacent slightly flattened hexagons, irrespective of the density, with the flattening due 
to the underlying square lattice. Note that, owing to the periodic boundary conditions, 
there will be configurations which appear to be different, but are equivalent to one another 
as they are connected by translations and/or rotations of the system. In the present work, 
such configurations are identified as equivalent and counted only once. 

Another issue that needs to be considered when placing the networks on a lattice is the 
mesh size dependence of the results, which should of course be negligible. Clearly, halving 
the lattice constant will lead to more states. But as long as the width of the energy interval 
SE used in the investigation is such that these new configurations all lie within the interval 
[E, E + SE], the change in the lattice parameter just adds a constant to the entropy. This 
corresponds to a parallel shift in a semi-logarithmic plot of the local density of states, and, 
as we shall see, does not affect our conclusions. We have ensured that this requirement 
is fulfilled reasonably well for our on lattice networks. Finally, one has to establish the 
connectivity of the configuration space. In our case, the neighbours of a given configuration 
are obtained by all possible moves of a single atom from its current position to one of the 
neighbour points of the lattice. With Na atoms we have in a d dimensional square lattice 
2dNA neighbors for each configuration. There is of course an amount of arbitrariness in the 
choice of the elementary moves which define the connectivity of the landscape. We were 
guided by the simple physical consideration that a single move should involve a change of 
coordinates which is small, i.e. of order a. This seems reasonable considering that collective 
moves are associated to vibrational motion, and thus take place on time scale much shorter 
than those of interest here. 



B. Lid algorithm 

The method chosen to investigate the energy landscape of the networks is the so-called 
lid algorithm. Here, we will only give a short summary of the basic procedure; a more 
detailed description of the algorithm and the problem-independent implementation we have 
used can be found in the literature |j,20|. Note that the lid algorithm is only applicable 



for discrete configuration spaces; continuous energy landscapes require a modified approach, 
the so-called threshold algorithm [I0|,p| . 



The central idea is to restrict the investigation of the configuration space to a set of 
smaller subregions, called pockets, which surround local energy minima. These pockets con- 
tain a few hundred to a few million states and can be explored exhaustively. The procedure 
is as follows: Starting from a minimum Xi, we list all the states that are accessible to a 
nearest neighbor random walk, which starts at x^ and which is restricted to states of energy 
lower than a prescribed energy value henceforth called the 'lid'. (By nearest neighbor 
random walk we mean a random walk in configuration space, whose steps consist of precisely 
the previously defined elementary moves). From the exhaustive listing, we can compute the 
the local densities of states D(E; Xj) and the local density of minima Dm{E; Xj) for 
states available below the lid L^. Integrating D(E\ x{) with respect to E over the inter- 
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val [0, Lk] we obtain the number of available states N(Lk]Xi), and by the same operation 
applied to Dm{E; Lk, Xj) the number of available local minima M(Lk]Xi). The search is 
repeated for successively higher lid values L k+ i, L k+ 2, etc., and for as many other pockets 
as possible. Note that the pocket is delimited by energy barriers instead of some distance in 
configuration space from the starting minimum. This is a natural first choice of selecting a 
physically relevant region of configuration space, due to the presence of an Arrhenius factor 
r esc oc exp(Lk/T), in the time of escape out of the pocket. One would like to repeat the 
analysis starting from every local minimum within a given pocket. Since the number of local 
minima within a given pocket of decent size (> 10 5 states) was already on the order of 10 4 , 
this was not possible, and the available computer time was instead used for a sampling over 
more pockets for a given cell size and density. The starting minima for these pockets were 
found by using either simulated annealing or an adaptation of the lid method itself. 

III. RESULTS 

In this section we first describe in some detail the features of the data obtained by the 
exhaustive investigations, and then discuss their links to the thermodynamical features of 
the model and of glasses. 

A. General features of the data 

Figures 2a, 2b and 2c describe three different pockets, belonging to the systems (Na = 
22; S = 12), (Na = 36; S = 16) and (N A = 53; S = 20) respectively. In each case we plot, 
on a semilogarithmic scale, N(L) and M(L) as functions of the lid L, and D(E; L max ) and 
Dm(E; L max ) as function of the energy E. For convenience, the dependence on L max of the 
last two functions is left understood. One notices the average exponential growth in all four 
quantities, with a flattening of D(E) and Dm{E) for energies near the maximum lid. Such 
a behaviour is exhibited in a large majority of the pockets, independent of (linear) size S 
and density p. 

Figure 3 shows semi-logarithmic plots of the local density of states D(E) for five different 
pockets belonging to systems of different cell sizes S and different densities p as indicated in 
the caption. These data are meant to illustrate the variations in the shape of D(E) seen in 
the simulations. Note that the larger systems shown in the panel a) have a steeper growth 
of the local density of states than the smaller systems displayed in panel b). In the following 
we denote by X{ the state of lowest energy within a certain pocket, and similarly by the index 
% various other quantities associated with the pocket. In about 60% of all pockets studied 
the curves for D(E) resembled those shown in (al) (26%) or (a2) (34%), i.e. they exhibited 
simple exponential growth D(E) = g i exp(E/E i ) = p i exp(a i £'), with g { = 1 and g; t > 1, 
respectively. Here, Ei = l/atj is the energy scale characterizing the exponential growth 
laws. In addition, the closely related case (bl) appeared for 29% of the pockets, while (a3) 
and (b2) represent some less common examples that occured in 3% and 8% of all cases, 
respectively. In some instances ((bl) and (a3)), the curves are best described by splitting 
them into two sections, each showing exponential growth, but with different 'growth factors' 
Q!j(l) and cxi(2). 
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Figure 4 is a plot of a =< Ei > _1 versus the linear size of the cell, for different densities. 
Here the averaging is performed over all the pockets of systems with the same density p 
and cell size S. It is clear that a increases and correspondingly < Ei > decreases with 
decreasing density and increasing cell size. This result agrees qualitatively and to a certain 
extent also quantitatively with what one would expect from the simple free volume analysis 



T4 outlined in the appendix. To show the agreement, the calculated curves based on eq.[ 



for p = 0.13, 0.14 and 0.15 are also depicted in figure 4 together with the actual data. 

For a given value of and S (considering only the larger pockets), one finds a consid- 
erable spread in the values of Ei, which vary up to a factor of two. This is quite different 
from the corresponding results obtained e.g. for spin glasses pip^fl , but it is consistent 
with the fluctuations of N(L) and D(E;L) around the average exponential growth curve. 
The size of the pockets and their local DOS are quite variable, possibly because different 
spatially localized excitation patterns of the network can exhibit different growth behaviour 
in the associated DOS. Thus large side-basins with different growth laws can appear upon 
exceeding some energy barrier. 

The density of minima shows a strong similarity to the density of states. Again, expo- 
nential growth is found, with the ratio of the growth factors cti(D) / 'a.i(D m) mostly in the 
range 1-2. This result agrees with the observation that the number of accessible states in 
a pocket N(L) is more or less proportional to the number of minima M(L), with the ratio 
N(L m ax)/M(L m ax) mostly in the range 15 — 60. 

As a function of L, N(L)/M(L) increases only slightly with the lid value. This holds true 
for all pockets investigated; and this observation is reflected in M(L) running nearly parallel 
to N(L) in a semi-logarithmic plot, once a lid value is reached where the first side-minima 
are accessible. Of course, D M (E; L max ) eventually decreases, in most instances for energies 
close to the lid, E ~ L max . 

It is natural to investigate to what extent the properties of a pocket depend on the energy 
of its lowest energy state. Denoting as the depth of a pocket the smallest energy barrier 
which must be crossed in order to gain access to a lower minimum, we find that at least 
for larger pockets, (N(L max ) > 1000), the depth grows as the energy of the lowest state 
decreases. With regard to the various growth factors, there exists a weak trend, insofar as 
pockets around low- energy minima tend to grow on the average slightly faster than the 
high-lying pockets, i.e. they have lower values of the E^s. However, many smaller pockets, 
often containing less than a dozen states and only a couple of minima, are interspersed 
among the larger ones along the energy-axis. Their depths do not appear to follow any 
particular pattern, but their growth quite closely follows an exponential law, analogous to 
the larger pockets. 



B. More detailed description of a typical example 

Let us consider in some more detail a typical case such as the system with (Na = 39, S = 
17). There are many local minima, each identified by its depth and by the energy of its 
lowest state. Each valley has its local density of states, which only includes states accessible 
by crossing energy barriers lower than the depth. All these local densities of states are to 
a good approximation exponential, and thus characterized by a growth rate, and its inverse 
Ei. Figure 5 shows a set of E^s (circles) and of depths (squares) plotted as function of the 
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energy E of the lowest state in the corresponding pockets. We see that the inverse growth 
rates do not vary much from pocket to pocket. The average value of E^ is for this system 
E™ ~ 0.035 ± 0.005 eV/atom. On the other hand, we note that pockets surrounding states 
of very low energy tend to be deeper than those surrounding less deep minima. In other 
words, the lower the energy the more rugged the landscape appears to be. 

In addition, we show in fig. 6 for the "ground state"- pocket ( "ground state" = deepest 
minimum found on the energy landscape) the differences 5DOS(E; L&) between the densities 
of states for subsequent lid values L k _i and L k . We note that these densities of states 
belonging to the added sub-regions at lid L are very similar up to the point of joining the 
main pocket. In particular, we note that at the lid L = 0.14 eV/atom a large group of 
similar basins that are about as deep as the starting minimum join the main pocket. If one 
considers the growth factors of these SDOS, one finds E^ v {5DOS) ~ 0.027±0.03 eV/atom. 
These values are quite similar to the growth factor of the whole "ground state"- pocket 
(Ei « 0.028) eV/atom, and lie at the lower end of the range of growth factors shown in 
figure 5. Thus, the larger sub-basins added with increasing lid-size are quite similar to the 
pockets encountered during the sampling of the whole energy landscape. 



C. Configurational entropy 

As mentioned in the introductory section, it is reasonable to discuss the excess entropy of 
the glass in terms of configuration space properties of network models. We have to assume 
that, at the temperatures of interest, the model system would thermalize in the pockets 
described by our numerical investigation. In this context, it is important to realize that the 
system experiences a qualitative change of behaviour at the temperature T — E^ Analyzing 
the expectation value of the energy of a pocket with an exponentially growing density of 
states D(E) oc exp(E/Ei), one finds that for T < Ei, the system is trapped in the local 
minimum at the bottom of the pocket, i.e. the high barriers of the pocket keep the system 
isolated from the rest of the energy landscape. But for T > E iy the system leaves the 
pocket with overwhelming probability, irrespective of the depth of the pocket. For a further 
discussion of exponential trapping and the competition among several exponential traps, see 



ref. [|3| 



If the assumption of thermalization within the pocket holds true for T < Ei, all quasi- 
equilibrium properties of the system can be calculated by the usual formulae of statistical 
mechanics, but with the sums over states restricted to a pocket. As the local density of 
states is available, we can calculate average energies and heat capacities. The model specific 
heat was calculated for the example pockets shown in figure 3, and plotted as a function of 
T in fig. 7. As one would expect, Cy shows a clear maximum at a temperature close to the 
average inverse growth factor of the local density of states. The height of this peak is the 
larger, the larger the pocket is, and it is more pronounced the less the DOS deviates from an 
ideal exponential growth law. Note that deviations from a perfect exponential growth show 
up as additional features in Cy. In particular, a rapid exponential increase (with Ei(l)) of 
the DOS at low energies followed by a slower exponential growth (with Ei(2) > Ei(l)) (c.f. 
curves (bl) and (b2) in fig. 3b)) is reflected in a prepeak at T « Ei(l) followed by the major 
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peak at a temperature somewhat below Ei(2) (curves (c) and (e) in fig. 7, respectively). . [] 
A qualitatively similar behaviour of the excess specific heat C p is observed in a large 
number of glass-forming systems, ranging from molecular and polymer systems to metallic 
and covalent glasses |f[4}| . Thus, the exponential growth of the local DOS of the network could 
be responsible for the configurational entropy observed in experiment. As a consequence, 



the glass transition would be the result of the system experiencing exponential trapping [23 



with the glass transition temperature Tq ~ Ei. In this context, one should point out that 
the high- temperature tail (T > Ej) of the calculated specific heat of a pocket has no direct 
physical relevance when comparing the model with real glasses, since for T > Ei the system 
would rapidly leave the pocket, and its equilibrium properties would no longer be dominated 
by the local density of states belonging to a single pocket. 

This hypothesis of the thermodynamics of the glass transition being controlled by the 
trapping temperatures of locally ergodic, exponentially growing regions of the energy land- 
scape of the glass raises the important question of the behaviour of the model system in 
the thermodynamic limit, N^, V — > oo with p = N^/V = const. Clearly, the number of 
possible neighbours of a configuration grows to infinity, and similarly the growth factor of 
the local DOS, i.e., the trapping temperature Ei goes to zero. This follows from the fact 
that due to the short range of the covalent interactions the energetic barriers in the system 
would be expected to grow only with y^ -1 )/^ while the energy, as an extensive quantity, 
grows proportionally to V itself. But one must not overlook the fact that there will be large 
entropic barriers preventing the system from exploring this infinite set of neighbouring states 
even though these are not separated by unsurmountable energetic barriers. From a certain 
point on, the system size has grown to such enormous proportions that the dynamics is 
controlled by entropic barriers, i.e. one is no longer allowed to assume that on the relatively 
short time scales available for observation the system can e.g. "focus" all the thermal energy 
present in the network into a precise sequence of moves needed e.g. to cross some barrier 
to a neighbouring basin. Else, we would be on time scales where the glassy state can be 
transformed into the crystalline one, with the consequence that the configurational entropy 
vanishes, of course. 

The existence of entropic barriers has important consequences for the dynamical be- 
haviour of the system. Since the trapping temperature is a local equilibrium quantity of 
an exponentially growing region 1Z of the energy landscape, it is necessary, in principle, to 
establish local ergodicity [0,0 within TZ, at temperatures below the trapping temperature. 
Usually, energetic barriers serve to delimit such regions, but in the thermodynamic limit en- 
tropic barriers fulfill this task. However, visualisation and characterization of such entropic 
barriers is usually not straightforward. 



1 This result agrees with the analytical calculation [^3| for the specific heat of a system restricted 
to a pocket, with an exponential DOS of depth Di between the minimum and the top of the 
(exponential part of the) pocket and an inverse growth factor E% : Cy = 1 for T << Ei, Cy = 1/T 2 
for T » Ei, and Cy = Di/12Ef for T = JSj. Note that the behaviour of Cy for T » Ei is a 
consequence of considering only the states with energies below the maximal lid L max . If e.g. the 
exponential growth is followed at higher energies E > L max by a power law growth, Cy approaches 
a constant value with increasing T after peaking at T ~ Ei. 
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The simplest picture of entropic barriers with regard to the energy landscape would be 
to assume that the configuration space of the excitations of the (infinite) system can be ap- 
proximately separated into a direct product of independent subspaces. Each such subspace 
can be treated in analogy to the independent modes of e.g. vibrations, and would usually 
be visualized as a "cluster of atoms" within the network that has its individual excitation 
spectrum and corresponding density of states. Such clusters would be only weakly correlated 
with each other. Since the actual number of states within a pocket of the landscape of such 
a cluster is small compared to the number of clusters in an infinite system, the distribution 
of energy throughout the system equals a Poisson distribution over the clusters. In partic- 
ular, the entropy S(E, V, N) becomes an extensive function of the number of clusters Nq'- 
S(E, V, N) = N C S C} where E c = E/N c , Vc = V/N c , and N AC = N A /N C ). The quantity 
Sc is the entropy of a single cluster, and for an exponential density of states is proportional 
to Ec/Ei, similarly to the examples presented in this paper. If one now assumes that the 
glass transition is a consequence of the energy landscape consisting of (possibly nested) 
locally ergodic pockets (e.g. the "clusters" discussed above) with exponentially growing 
densities of states, one would identify Tc with the average inverse growth factor Ei(Vc) of 
such regions. Using the simple growth law (Ei oc 1/V) derived in the appendix, such an 
identification yields an estimate of the size of these clusters. Since in the independent cluster 
approximation clusters of size Vc suffice to describe many of the thermodynamic properties 
of the networks, this establishes a reasonable network size for structural investigations at 
the level of network topology. 

In the case of 2d-networks, the trapping temperature for a 20 x 20 system with 50 — 60 
atoms lies in the range of about leV(« 10 4 K). Since the glass temperature for covalent net- 
works lies in the range of 0.1eV A (~ lO 3 ^), the energy landscapes of networks with N A > 500 
might be suitable for a realistic description of properties of glasses. While this number has 
been derived from the analysis of 2d-networks, preliminary results for 3d-networks indicate 
that the trapping temperature for e.g. networks with 40 atoms lie at about 1.2 eV/atom, 
leading again to an estimate of the needed network size of about Na > 500. 



IV. CONCLUSIONS 

In this paper, we have presented a first in-depth study of the microscopic energy land- 
scape of glasses at the level of their network topology. Such amorphous networks show an 
exponential growth of many important quantities, in particular the local density of states 
and minima, D(E;L) and D M (E;L) respectively, the accessible state space volume N(L) 
and the accessible neighbor- minima M(L). This growth leads to exponential trapping that 
can explain the occurrence of a glass transition in such systems. In particular, the peak in 
the excess specific heat at the glass transition, that is often observed in experiments, follows 
directly from the exponential growth law. 

This type of local exponential growth appears to be a common feature of many complex 
systems. Similar behaviour has been found for polymers flr2" , |55 |, spin glasses pl| , p2| , |26| , 



combinatorial optimisation problems || , crystalline solids p7| , |28| and, in preliminary results 
by the present authors, for three-dimensional random networks. 

It is expected that these results can form a more solid basis for phenomenological models 
of complex systems. An important open question is the issue of entropic barriers, since they 
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control the effective size of locally ergodic pockets in large systems. Knowing this size would 
allow the direct quantitive comparison of the calculations with experimental data; e.g. one 
could predict the glass transition temperature and its dependence on cooling rates. 
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V. APPENDIX 

For the qualitative and semi-quantitative analysis of the network glasses, the following 
very simple free- volume- mo del that is based on some rough assumptions about the nature 
of the energy landscape can be helpful: 1. Starting from each configuration, 2dN A new 
configurations can be generated by moving each atom to one of the 2d neighboring sites on 
the lattice (d is the dimension of the lattice). 2. The interactions are highly local, such 
that the typical increase in energy Em of an acceptable move, i.e. one which leads to a 
configuration below the energy lid, is essentially independent of the number of atoms in the 
simulation cell. The typical increase in the energy/atom for a system with Na atoms when 
an acceptable move is performed is thus E m /Na- 3. As long as the states belong to the fast- 
growing region of the "pocket", the number of downhill moves among the acceptable moves 
will be very small. Thus we assume that a given move will result in a new configuration with 
an energy increase SE ~ Em with a probability f(p, Na)- In particular, the probability /_ 
for a downhill move should be very small compared to /, /_ << /. Thus, the number of 
configurations N(E + 5E) with an energy below E + 5E, is proportional to the number of 
configurations below the energy E, N(E): 

N(E + 5E) oc N(E) (2dN A f) (3) 

from which it easily follows that 

N(E) = N exp(aE), (4) 

with a = l/Ei = (1/ ' E M ){2dN A f — 1) ~ 2dN A f / E M - The last approximation should be 
acceptable in the limit Na —> oo , as long as / varies only weakly with N A (for fixed p). 

It seems reasonable to assume that / depends on the density of the system, and that it 
increases with the free volume per atom Vf = Vf /Na in the system, where the free volume 
Vf is given in terms of the volume per atom Va as Vf = V — N a Va- Thus we get, for some 
constant c 

f = cv f = c(V - N A V A )/N A = c{V/N A ){\ - pV A ), (5) 

and therefore 

a = l/Ei » 2dN A f/E M oc V(l - pV A ) (6) 

Equation |^ shows that / increases with decreasing density for constant volume. Thus 
also decreases with density, and secondly, for constant density 
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Ei oc l/V. 



(7) 



This final result agrees quite well with the observed behaviour. 

Note that according to eq.|, Ei should actually decrease as 1/V(1 — pVa). In order to 
check this behaviour, we need to determine Va- This value corresponds to the number of 
lattice points that the atom would occupy in an energetically favourable configuration (recall 
that we are investigating the pockets around deep- lying minima). For an approximately 
hexagonal arrangement on the square lattice, we find Va = 6. If one plots Y{= NaE^) 
against X{— 1/(1 — pV^)) for all available volumes V = S 2 (S = 10, ...,20) one finds that 
the data roughly follow straight lines as suggested by eq.|6| (T2| . 

Note that there exists only one free parameter (c/Em) in this simple free volume model. 
The number of accessible states N(L) estimated by this model exhibits an average expo- 
nential growth similar to D(E) (c.f. fig. 2), and the qualitative dependence of the growth 
factors on V = S 2 and p should be essentially the same for D(E) and N(L), unless e.g. 
there is a high degree of degeneracy in the ground state of the pocket. But even for polymer 
systems, where such a degenerate ground state occurs more frequently, an analogous free 
volume model roughly describes the dependence of the growth factors on system size IE 
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FIGURES 



FIG. 1. A typical network configuration belonging to a low-energy metastable state in a 
system of 27 atoms with a unit cell of linear size 14. 
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FIG. 2. Number of accessible states N(L) (squares), accessible minima M(L) (circles), density 
of states D(E) (triangles) and density of minima Dm{E) ( diamonds) within a pocket, as a function 
of L [eV/atom] and E [eV/atom], respectively. Data are shown for three different pockets: a) 
(N A = 22; S = 12), b) (N A = 36; S = 16) and c) (N A = 53; S = 20). 
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FIG. 3. Local density of states D(E) vs. E [eV/atom] for several representative systems: a): 
black squares (al): (Na = 46; S = 18), black circles (a2): (Na = 41; S = 17), and black triangles 
(a3): (N A = 34; S = 16). b): white squares (bl): (N A = 13; S = 10) and white circles (b2): 
(iV^ = 21; 5 = 12). Note the curves belonging to the larger systems shown in a) are much steeper 
than those belonging to the smaller system shown in b). 



15 




10 11 12 13 14 15 16 17 18 19 20 

linear cell size 



FIG. 4. The quantity a = 1/ < Ei > [l/(eV/atom)] for different densities: black squares: 
p = 0.13, black triangles: p = 0.14 and black circles: p = 0.15. The white symbols show the 
corresponding theoretical predictions obtained through eq|| with Va = 6. The average over the 
inverse growth factors Ei is performed over all pockets in systems of the same density and cell size. 
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FIG. 5. Depth (squares) and inverse growth factors Ei, both quantities expressed in [eV/atom] 
for several sample pockets as a function of the lowest energy of the pocket. The system considered 
has (N A = 39;S = 17). 
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FIG. 6. Local densities of states of several sub-basins of the 'ground state' pocket in the 
system (Na = 39; S = 17), which appear successively as the lid is raised. The top curve is the local 
density of states of the whole pocket. Note that the exponential trend is common to all sub-basins. 
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FIG. 7. Specific heat in units of [ks] vs. temperature [eV/atom] calculated for several rep- 
resentative examples: a) (N A = 46; S = 18), b) (N A = 41; 5 = 17), c) (N A = 13; S = 10), d) 
(Na = 34; S = 16), and e) (Na = 21; S = 12). The three peaks in plots a), b) and d) are located at 
temperatures which match the inverse growth factors E\ of the exponential density of states shown 
in Fig. 3a), (al), (a2) and (a3), respectively. The double peaks of plots (c) and (e) correspond 
to the two different slopes characterizing the densities of states shown in Fig. 3b), (bl) and (b2), 
respectively. 
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